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1 Introduction 



The formation of structure in the Universe is one of the most interesting mysteries 
in nature. The latest observations Q have revealed that for galaxy formation, density 
perturbations must grow much faster than the Universe expands. However, the standard 
theory of gravitational instability in Einstein's gravity theory does not provide such a 
high growth rate. Although some modifications of the standard theory with appropriate 
dark matter, a cosmological constant, or cosmic strings have been proposed, we do not 
have any completely satisfactory solution yet. Hence it is worthwhile to search for an 
alternative, in which the gravity sector is modified. 

Generalizations of Einstein's theory have found various interesting applications to 
astrophysics 0-0. The Brans-Dicke (BD) theoryQ or its modified versions may provide 
a graceful exit from the inflationary era of the Universe (Extended or Hyper-extended 
inflation [Q], and soft inflation^]). Non- minimally coupled terms may produce a virtual 
structure in the Universe ||. There may also exist an important difference in the growth 
rate of cosmological density perturbations, which might resolve the problem of structure 
formation. In a preliminary work, Berkin and one of the present authors investigated 
the growth of density perturbations in a theory with a non-minimally coupled scalar 
field and showed a sufficient enhancement in the growth of scalar field perturbations^. 
Although these authors excluded ordinary matter fluids and used only a scalar field, 
their calculation suggests that such a coupling may play an important role in structure 
formation in a realistic model with a baryonic matter fluid. 

In this paper, we discuss cosmological density perturbations in generalized Einstein 
theories (GETs), which include the BD theory, theories with a non-minimally coupled 
scalar field and certain curvature-squared theories. The action of GETs we consider here 
is 



S = d 4 x\f^g 



2 

m 



XI) 



where F((p, R) is an arbitrary function of a scalar field (f) and a scalar curvature R, e(4>) 
is an arbitrary function of (f), and L m is the Lagrangian of ordinary matter, for which we 
assume a perfect fluid. 

When we study cosmological perturbations, an appropriate choice of gauge becomes 
important. In order not to pick up unphysical modes by a bad gauge choice, Bardeen pro- 
posed a gauge-invariant (GI) formalism for density perturbations. Although Bardeen's 
GI formalism H is one of the most elegant and attractive approaches to cosmological 
density perturbations, its application to GETs requires tedious calculations, which have 
yet to be performed]]] 

On the other hand, Ellis and his collaborators have recently developed an alterna- 
tive to Bardeen's formalism [pTf - |l4|[ . Their covariant approach is rather simple. Using 
their formalism and a conformal transformation [JD]]- [16], we will present GI perturbation 

1 Kodama and Sasaki nicely reviewed and extended Bardeen's formalism in p0| . They also 
tried to formulate GI equations for a theory with a non-minimally coupled scalar field, but owing 
to the complexity of the problem, they made no attempt to solve it completely. 
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equations in GETs. Although we can write down the basic equations schematically for 
an arbitrary GET (see appendix A), we will present them explicitly only for the following 
two interesting cases: 

Case (A) : F is a linear function of R, i.e., F{<j>, R) = f((f>)R - V{4>) (in §3). 

Case (B) : F depends only on R, i.e., F((p,R) = L(R) (in appendix A). 

Case (A) includes the BD theory, induced gravity[i~7|l, and other theories with a non- 
minimally coupled scalar field, while Case (B) includes -R 2 -gravity||]. 

In order to show the possibility of an enhancement in the growth rate of density 
perturbations in GETs, we analyze perturbations in the BD theory. 

The plan of this paper is as follows: A brief review of the covariant approach to den- 
sity perturbations is given in §|2[ In §|| we derive the perturbation equations for Case (A), 
introducing new conformal variables. In §|| we apply these equations to BD cosmology 
and analyze the evolution of density perturbations of dust fluid both analytically and 
numerically. Finally in §|| we discuss general features of density perturbations in GETs. 
Appendix A presents the formulation for the most general type of GETs described by 
equation (1.1) and gives explicit expressions for Case (B). In appendix B we show our 
equations reduce to the already-known ones in a scalar field dominated universe. In 
appendix C perturbations in a BD cosmology are solved analytically. 



2 Covariant Approach to Cosmological Density Perturba- 
tions 

The covariant approach to cosmological density perturbations, introduced originally by 
Hawking |]18|| and developed by Ellis and BruniQ], provides us with a GI formalism, 
which is an alternative to Bardeen's approach p|-|10||. Here we summarize it briefly in 
order to define our notation and introduce the dynamical variables used in this paper. 

In the formalism, first we have to select a frame in which GI quantities are defined. 
The frame is described by 4- velocity of an observer, u a . The time derivative along the 
motion of the observer (denoted by a dot) of any rank of tensor T ab ' cc i... and its spatial 
covariant derivative in the 3-space S orthogonal to u a (denoted by ^ 3 ^V a ) are defined by 

t" h ,/V, •/•"'' ,,/..- (2.1) 

and 

::i T„T>" = h a m K h hf ■ ■ ■ h/hj ■ ■ ■ V,„V l (2.2) 

respectively, where h a b = g a b + u a Ub is the projection tensor into the 3-space X and V a 
denotes the covariant derivative with respect to g a b- 
The derivative of u a is decomposed as 

VbU a = Wat + <?ab + \oh ab - a a u b , (2.3) 

where a a = u a , 6 = V a ii a , uj ab =^ V[ b n a j , and a ab =^ V( fc u a ) are the acceleration, 
the expansion, the vorticity tensor and the shear tensor, respectively. We define typical 



2 



length scale a and time scale H from the expansion 8 as 

\e=± = H. (2.4) 
6 a 

These a and H turn out to be the scale factor and the Hubble parameter, respectively, 
when the Universe is homogeneous and isotropic. In inhomogeneous spacetimes, however, 
those are local quantities defined by each observer. 

Next we define the GI perturbation variables. The basic requirement for GI quan- 
tities is that they are invariant under a general coordinate transformation, i.e., for any 
choice of correspondence between a homogeneous and isotropic background spacetime 
and the physical, inhomogeneous universe. The simplest GI quantities are a scalar field 
which is homogeneous in the background, and any vector or tensor field which vanishes 
in the background. For such quantities, all perturbations defined by comparison with 
background values are free from any gauge dependence. Such quantities as a a , u a b and 
a a b defined in a real lumpy universe are GI perturbations, because their background 
values vanish. 

As for the energy density fi, the spatial variation 

X a = ( 3 )V a /i (2.5) 

is GI because fi is homogeneous in the background. The GI quantity 

V a = -X a , (2.6) 
I* 

which is the ratio of its spatial gradient to the density at a fixed comoving scale, is con- 
venient in discussing cosmological perturbations. If we are interested simply in density 
perturbations, however, we have to extract the information of its scalar part from the 
vector quantity T> a , which contains other information as well. In this formalism, a local 
and unique splitting Jl^] is attained by the operation of the (comoving) spatial derivative 
a^Vj, and its decomposition. For T> a , we have 

a^V b V a = A ab = ^Ah ab + Z ab + W ab , (2.7) 

A = h ah A a b W a b = A[ ab ], Ti ab = A( a6 ) — -Ah ab , 

where the skew-symmetric tensor W ab contains information about vorticity, the symmet- 
ric and trace-free tensor T, ab describes the evolution of anisotropy in the universe, and 
the trace A is related to the aggregation of matter fluid. This A, which corresponds to 
the density perturbations Sfi/fi, is one of the most important GI variables^. 

2 A is related to Bardeen's GI variable e m (= Sfi/fi in velocity- orthogonal slicing) in || as 
follows @: 

n n 

to first order of perturbations, where n is a wave number of the harmonics Q^ n \ 
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Obtaining the evolution equations for GI variables is rather simple. First, take 
spatial derivatives (by operating a^V a ) of the fundamental equations, which consist 
of the energy conservation, the momentum conservation, the Raychaudhuri, and the 



Gauss-Codacci equations (see e.g. [19] or ]p0|| ). Then, take the (comoving) divergence 
(by operating a^X7 a to the gradient equations) and construct scalar equations for GI 
variables such as A defined above. The resultant equations describe the evolution of GI 
scalar perturbations. 

3 Formulation of Perturbation Equations in GETs 
3.1 Conformal Transformation 

We now derive GI perturbation equations in GETs, using the covariant formalism de- 
scribed in §^. In this section, we consider Case (A), i.e., the case where F(<p,R) in 
equation ( |1.1|) is a linear function of R. As for Case (B), i.e., the case where F(c/),R) 
depends only on R, we will present explicit perturbation equations in appendix A, which 
turn out to be the same as those given here. Furthermore, the derivation of perturbation 
equations for the most general action (1.1) is straightforward (see also appendix A). 



The theories we consider in this section are described by the action 



S = / d A x^=g~ 



/(</»)i?-^(v«/») 2 -y(0) + L r 



(3-1) 



where f(4>), c(<j>) and V{4>) are arbitrary functions of cf>. L m is the Lagrangian of ordinary 
matter, for which we assume a perfect fluid, i.e., 

m = l*v>fvF+ph£, (3.2) 

where /z, p and u^f represent the energy density, the pressure and the 4- velocity of the 
matter fluid, respectively, and is the projection tensor into the 3-space orthogonal 
to u*f . Taking variations of equation (|3.1|) with respect to the metric g a b and 4>, we find 
the basic equations 



Gab 



+ [v a v 6 /(0) - n/(0) ff , 



-\v{4>)gai> + \T^Y (3.3) 



3D/(«£) + ^fW) 2 + 2V(</>) + - 3p)} = 0, (3.4) 



where G a b = R a b — \§abR is the Einstein tensor. 

The application of the covariant approach in §2 to the present model is not straightfor- 
ward because of the presence of higher-derivative terms such as V a V&/(</>) and Cf(4>)g a b- 
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Although there may be a way to resolve this problem leaving the present variables as they 
are, we prefer to circumvent it by changing variables using the conformal transformation (l? 

9ab = e 2 ^g ah . (3.5) 

If we choose the conformal factor as 

e^ = 2K 2 \f{4>)\, (3.6) 

the field equations become 

ee K 2 (f*+f%), (3.7) 
□ <p - U'((p) + n 2 M(ip){n - 3p} = 0, (3.8) 



with 



e(^)/(0)+3(d/(^)/#)- 
2/ 2 (0) 



(3.9) 



u ™ = WJW V ^ (sign) = ^W' (3 - 10) 

M(cp) = r> (3.11) 

(2 K 2/(0)) 2 [2K 2 {e(0)/(^) + 3(df (0)/#) 2 }]^ 

N M = wfm> (3 - 12) 

= e"uf, u* M = e~"u* M , (3.13) 
= K\ h% = e 2 "h%, ht = e~^h% (3.14) 

where k 2 = 8irG, ' (a prime) =d/d(p, and new variables with A (a caret) denote those 
with respect to g a b- The resultant basic equations in terms of new variables turn out 
to be familiar ones, that is, those for the Einstein gravity (gab) plus a minimal scalar 
field tp with a potential U((p) and a perfect fluid, which interact with each other through 
coupling functions M((p) and N((p). The non-minimal coupling to a scalar curvature has 
been absorbed into such interactions. Note that we have not changed the original energy 
density fj, and pressure p for a perfect fluid. 

We can recover physical variables from new ones through equation ( [3.5D , which gives 
the relations between the length scales and Hubble parameters 

a = e w o, (3.15) 
H = e~ u (H + u), (3.16) 
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and equations ( |3.14 ). In particular, the density perturbation A M observed along fluid 



flow u„ is obtained from A as follows: First, we have the relation 

V M _ °(3) V M = £ fc Af6 v = e -o, £^Af6^ = e - W «(3)^ _ e -u^M (3 . 17) 
/' /' /' /' 

to first order in the perturbations. Taking the comoving divergence of this relation by 
a^V a M (= e u a^V a M ) yields 

A M = A M , (3.18) 

Thus, there is an advantage to using the conformal variables. The original covariant 
approach is applied as it is to the simplified conformal equations and it gives the same 
density perturbations as those in terms of the physical variables. In what follows, we 
will formulate equations in terms of the new conformal variables. We omit * (a caret) 
for brevity. When we discuss the results, of course, we transform back to the physical 
variables (see §4). 

3.2 Formulation in terms of Conformal Variables 

The total energy- momentum tensor, T* b = T^ b + T^ b in equation (|3.8|) , takes on a fluid 
form for any observer with 4- velocity v,® [13], 



T*b = V*u a u b + p*h ab + q <a u b) + vr^, (3.19) 

with 

M* = ^ 2 + U(ip)+N(<p)n, (3.20) 

P* = \v 2 -U(cp)+N(cp)p, (3.21) 
q° a = -^V a ^ + N(^+p)V a M , (3.22) 
v?ab = 0, (3.23) 

where we have neglected the higher than first order of perturbations composed of all the 
spatial derivatives and the relative velocity of the matter fluid with respect to u^, 

V a M = uM-u° a . (3.24) 

In what follows, we will also neglect higher-order perturbations in all equations, as in the 
usual first order perturbation theory. In equation ( 3.19| ), /i*, p*, q® a and ir® ab represent 



the total energy density, pressure, energy flux and anisotropic pressure measured by the 
observer ■ Note that the energy density ^* and the pressure p* do not depend on the 
observer, but only the energy flux does. 

There is no physical constraint on u® other than that V^ 1 be small. However, there 
exists a preferred choice, i.e., the 4- velocity of the center of mass. In this reference frame, 
the total energy flux vanishes. Ellis et al. call it the energy frame and we use to 
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denote its 4-velocity. The total energy-momentum tensor, written in terms of the center 
of mass observer (u^), takes the perfect fluid form: 



T ab = fi*u^u b +p*ha b , (3.25) 

where //* and are defined by equations ( |3,20| ) and ( 0.21 ) . 

Once a reference frame is fixed, we can derive explicitly the evolution equations from 
energy conservation, momentum conservation, the Raychaudhuri equation, the equation 
of motion for <p, and the Gauss-Codacci equation in the center of mass frame, which are 

fi* + 3H(n*+p*) = 0, (3.26) 
(^+P*)a a + (3) VaP* = 0, (3.27) 
3H + 3H 2 - V a a a + -k 2 (h* + 3p*) + 2(a 2 - uj 2 ) = 0, (3.28) 

^R = 2(-\e 2 + ^ + a 2 -uj 2 ), (3.29) 

U(p-U'(cp) + K 2 M((p)(fj,-3p) = 0, (3.30) 

where a 2 = \o ab a a b and uj 2 = \uj ab uj a b. 

The equations for the background isotropic and homogeneous fields are derived by 
setting cr, u and the spatial derivatives equal to zero in equations ( |3.26| )^( 3.29 ) as 

A + 3H(n + P ) + + K '^y 0* ~ W = °> ( 3 - 31 ) 

3H + 3H 2 + K 2 (ip 2 - U(ip)) + ^K 2 N(<f)(fi + 3p) = 0, (3.32) 

3k ( 1 \ 
3H 2 + = k 2 i^-tp 2 + [%) + iV(^) M J , (3.33) 

<p + 3#<^ + - K 2 M((p)(n - 3p) = 0, (3.34) 

where k{= or ±1) is the so-called spatial curvature constant in the Friedmann-Lemaitre- 
Robertson-Walker(FLRW) universe. 

The scalar GI perturbation variables for the matter and scalar field in the center of 
mass frame are defined as 

Af = -( 3 >Vi^, y? = a 2 ^vW (3.35) 

A E = ^!(3) V | jU; yE = fl 2(3) v |, pj (3.3,5) 

^£ = a 2(3) V|# (3.37) 
^ = a 2{3) V 2 EV . (3.38) 
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Following the procedure in §2, we find the first order perturbation equations as 
follows: 



A 



o, 



4 + (3) V| )3f = 0. 



i$> E + 3^6 £ + {[/"(^)- K 2 (M'((^))(^-3p)}$ E 



(3.39) 
(3.40) 



- K 2 M(<p)^A E + 3 K 2 M(v)y E + - ( 3 ) V 2 E <S> E _ aV (3) V|a a 
-a 2 (4^ + 2(^) (3) V|o a = 0. (3.41) 

Here we have used the following useful relations for uj^, <p, an arbitrary vector T E 
orthogonal to u% (T E u% = 0) and T E = a^V%T E : 



a (3) V |(3) V ^ afc = 0, 
Uip = \$ E - 3Hip - tp, 



a^V a E T E = t E , 



a (3) v a fE = f £ 



2k, 



( (3) V| - ^)Tf 



(3) V | T S 5 



(3.42) 
(3.43) 

(3.44) 



to first order of perturbations. See the appendix of [|T^] for more details. 

The explicit equations in terms of A E , y E , & E and Z E are derived from equations 
(U|)~(^||), by using relations (|3.20|) , and 



A; 



l<p + U(<p)+N(<p)/A 



N(ip)^A E - a 2 ^V a E a a 



+^ E + (U'(<p) +nN'(<p)) & 



(3.45) 



Equation ( |3 . 2 7 ) is also used to eliminate the acceleration a a . To describe the resulting 
equations, we expand the dynamical variables as 

A = ^A (n) ( f )QW, (3.46) 

n 

where is an eigenfunction of the Laplace-Beltrami operator ( 3 )V 2 [14]: 



(3) v 2g(n) 



rr 



Q 



(n) 



(3.47) 



Finally, we obtain the basic equations for the harmonic components of perturbed 
quantities as follows (We drop the suffix (n) for brevity): 



A 1 



3w (h- k 2 ^) A e + 3* 2 ^y E - (1 + w) Z 1 
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+ 



+ 



3H 



(W - N'fiw) + p ■ 



~a 2 N~Ji 



~N 



M 



) (1-3^) 



(3.48) 



~K 2 NnA l 



1 



2 

1 

2 



K 2 + 



2fi(l + w 
1 / 



2 -2 

K (p 



2n2 

^2" 



6 A; 



z ^r + ^\y E - 2HZ E 



Nfi(l + w) 



2 -2 

K 



2n 2 



6fc 



a 2 a 2 



k 2 {U' + N'n) 



{V - N'nw) 
N/j,(l + w) 



2 -2 

K if 



3> e 

2n 2 
^ 2_ 



6A; 

+ ^7 



<S> E , (3.49) 



^ 2 + -/V>(l + U;) 



,M 



3{H + k 2 — p)N l i{l + w) + 



1 + c 



l + w 



V 



N' n M , 



N' 



—p 3 + p {3(c 2 - l)iT^ - W - 12k 2 M^w + N'fiw} 



ip 2 + iV/x(l + w) 



n 



U" - k 2 M'h(1 - 3u>) + ^ + 



(U" - N"fiw)p 2 + c 2 s ^-fip 2 {N' + k 2 M{1 - 3w)} 



3c 2 HN'p ' 



N 



Nfj,(l + w) 



-{U' -N'fiw)- 



N'w 



"s -2 



N 1 + w 



if — K 



,M 
~N 



(1 - 3w) ( 



2N/I + 



1 + W 



-(3c 2 s Hp-2U') 

K 2 MNft 2 {l + w) 
p 2 + Nfj,(l + w) 
1 

p 2 + iV/x(l + to) 



p 2 + Nfi(l + u; 



K 2 MNfi{h - 3w) + N(3c? s Hp - 2U' 



1 + c z s . 2 

1 + W 

Nfi(l + w)p 
p 2 + Nfi(l + to) 



{iV' + K 2 M(l -3w)} 



3^ 



(3.50) 



where we have used notations w = p/ fi, c 2 = p/fi^. 

These complicated equations may be much simplified in the case where N((p), M{ip) 
and/or U(ip) are expressed in terms of exponential functions of tp. This happens in 
several interesting models (see eqs. ( [4.15| )^(4.17) for the BD theory and also appendix 
A). 

Although the variables in the center of mass frame are easy to handle, when we discuss 
the density perturbations of matter fluid, we should evaluate the density fluctuations of 
the matter fluid in the matter frame defined by uff . Hence, we must construct the 



density perturbations in the matter frame from variables in the center of mass frame [13 
First, the projection tensor defined by uff and the energy density of matter fi M in 
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the matter frame are written as 



«a& - 5a6 + U a U b - tl ab + U a V b + V a U h , 



to first order. Then, using these relations and the condition 

V M = t 

aN{n+p) a ' 

which is derived from the condition that the total energy flux vanishes in the center of 
mass frame, the density perturbation in the matter frame A M is expressed as 

= A E + ±—^—<S> E . (3.51) 

In some specific theories, we recover already-known basic equations. For example, 
setting f{4>) = 1/2k 2 , and e{4>) = V{4>) = 0, we find immediately the perturbation 



equations for the dust universe in the Einstein theory. Hwang[21] and Bruni et al. [|22| 
derived perturbation equations in some scalar field dominated universes. We can also 
recover these equations by taking the limit /x, p — > and rewriting our variables in terms 
of their variables (see appendix B). 

4 Density Perturbations in the Brans-Dicke Theory 

In order to see how the coupling term to the curvature scalar in GETs changes the growth 
rate of density perturbations, we consider the BD theoryQ as an example. In the BD 
theory, the reciprocal of the scalar field, l/<p, corresponds to an "effective" gravitational 
constant G(4>). If a flat dust universe (k = 0, p = 0) satisfies a specific initial condition 

/Mo = (3 + 2a;) 

0o 4vr(4 + 3w)' 1 ' ; 

where a; is a constant (the BD parameter) and the suffix represents the present time, 
the BD scalar field 0, the energy density of matter /i and the scale a evolve in time as 



2/(4+3w) 

(4-2) 



to J 

h) ■ ,4 - 3) 

/M -3(2+2w)/(4+3w) 

H = n a- 3 = M i^-j . (4.4) 

In the limit of large u, this solution approaches the Einstein de-Sitter universe (4> = 
constant, a cx t 2//3 ). When the condition (4.1) is not satisfied initially, the behaviors of 
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and a in the early stage may differ considerably from the above solution for the same 



value of u>[23|. 



Nariai ||24| analyzed density perturbations in the above specific solution (4.2)~( fL4|) 
of the BD cosmology using a specific gauge condition and solved its evolution equations 
analytically as 



+ Q^j 2 {c 3 Ju(nT) + c 4 J-v(nT)}, 



(4.5) 



where c%, C2, C3 and C4 are arbitrary integration constants, v = (4 + 3a;)/2(2 + uS), 
T = 2vto(t I 'to)~ ^ , Jy is the Bessel function of the u-th. order and n represents the wave 
number of harmonics. 

In our formalism, the BD theory is given by setting 



0, 



v(<f>) 



(4.6) 
(4.7) 
(4.8) 



in the action (3.1). New variables, potential and coupling functions in the conformal 
frame introduced in §3 are given as 



dt 
df 



(4.9) 



N(ip) 

M{ip) 
[%) 



v/2(3 + 2lo) 



2k 



In 



exp 



4k 



^2(3 + 2u) 



k^2{3 + 2lo) 
0. 



Nip), 



(4.10) 

(4.11) 

(4.12) 

(4.13) 
(4.14) 



The exponential form of N and M and vanishing U rather simplify the perturbation 
equations ( 3.48 )~ ( [3.50 ). Moreover, if we assume the ordinary matter is dust (p = 0), 
they are reduced as 



+ 



3 k 



3Hip 

Nn ' ^2(3 + 2lo) 



n rip P 
a 2 N/i ' 



(4.15) 
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Z E 



k 2 N\iA e - 2HZ 1 



k 2 + 



2 -2 

K if 



2n 2 6/c N 



+ 



^2(3 + 2uj) 
1 



(4.16) 



+ 



ip 2 + Nfi 

N/j, 
ip 2 + Nfi 
kN 2 fi 2 



3H{Nfi - <p 2 ) + 



Kip 



V2(3 + 2w) 



(3iV/i + 



2k 2 N{i 
3 + 2lu 

1 



Vi 2 + AT/i ^/2(3 + 2a;) 



n 



A 1 



Nyip 
ip 2 + N/j,' 



(4.17) 



In the case of a flat dust universe (k = 0, p = 0) with the specific initial condition 
Q4.1D , we can solve these equations analytically (see appendix C). As for the density 
perturbation of dust, we obtain 



A 1 



+ c 2 




(4.18) 



where c\, C2, C3 and C4 are arbitrary integration constants, and 



4 + 3a; 
2(2 +lo) 



4 + 3uj 1 5 + 3a; t 



5 + 3ui 



4 + 3a; t 



(4.19) 



which are the same as Nariai's notation. We can easily show that this solution is essen- 
tially the same as the gauge-dependent one ( f4.5[ ) by use of the relation ( |C.ip , 

The density perturbation A E has four independent modes. The first two terms have 
their counterparts in the Einstein de-Sitter model in the limit of uj —* 00, while the last 
two originate from the BD scalar field. The growing mode (the first term in equation 
fl4.18| )) evolves asymptotically 



"'To 



2(2+u;) 
i+3uj 2+ui 



(4.20) 



which shows that the growth rate of A E in this specific BD model is somewhat higher 

1 

(by times) than that in the Einstein de-Sitter model. 

Although we have obtained an analytic solution for the above specific model, in 
order to see whether such an enhancement in the growth rate of density perturbations 
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is generic and to find what the key to such an enhancement is, we reanalyze numerically 
the density perturbations of BD cosmological models with arbitrary initial conditions 
using our formalism. 

We also evaluate the growth of density perturbations quantitatively under constraints 
from observations. Thus we can confirm that if all the observations are reliable, then 
the enhancement in the BD theory itself is inadequate to resolve the structure formation 
problem. The observational constraints on the BD theory may be summarized as follows: 

1) The BD parameter |25[| ; 

uj > 500, 

2) The present rate of time variation of (j|26[|; 



< 1 x 10~ n year" 



present 



3) From nucleosynthesis [ 27 1; 



AG 



G 



< 40%, 

nucleosynthesis 



where AG = G - G\ present . 

Among these constraints, 1) and 2) are obtained directly from present-day obser- 
vations of gravity, while 3) is a result obtained indirectly, i.e., through the comparison 
between the theoretically predicted and observationally inferred primordial abundances 
of light elements. If we adopt 3), which critically depends on the theory of nucleosyn- 
thesis, the time variation of gravitational constant in the past is also strongly restricted, 
resulting in no difference in the evolution of density perturbations from that in the 
Einstein theory, as we will see shortly. 

For the sake of easy comparison with the Einstein de-Sitter model, we first assume the 
ordinary matter is dust (p = 0) and the present density parameter CI dust, o = k 2 Ho/3Hq 
equals unity, which means the universe is closed (k > 0) in the BD cosmology. Besides, 
we set the Hubble constant Hq = lOOkm/sec/Mpc. Then, regarding the observational 
constraints, we consider the following two cases: 

Case (i) : All the constraints 1), 2) and 3) are satisfied. 

Case (ii) : Only 1) and 2) are satisfied. 
In all the calculations below, we set the BD parameter u = 500, therefore the constraint 
1) is always satisfied. As for constraint 3), we calculated the background spacetime back 
to the nucleosynthesis era where a ~ 10~ 9 and found \G/G\ present < 2.039 x 10~ 13 year 
for the universe which satisfies the constraint \AG/G\ nuc ieosynthesis < 40%. Hence we 
adopt \G / G\ presen t = 2.039 x 10 _13 year _1 for the model in Case (i). As for the model in 
Case (ii), we set \G/G\ present = 2.042 x 10~ 13 year _1 , resulting in | AG / G\ nudeosynthesis ~ 
200%. Both values for \G/G\ presen t clearly satisfy constraint 2). 

For the above two cases, we calculate the evolution of the density perturbations of 
dust A E as the scale of the universe grows 10 3 times to today. We assume the wave 
number of the perturbations n = 30, which means today's ratio of physical wavelength 
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of perturbations to the Hubble horizon scale Xo/Xh,o = 1/30, corresponding to large scale 
structure on scales of Ao — lOOMpc. We normalize the initial condition for A E (= Sfj,/fj,) 
to unity at the starting time. Therefore, it is natural to set the initial conditions for 
Z E (= 35H) and $> E (= dip) so as to approximate roughly their background values H 
and ip, respectively. We set Z E ~ -H(= -32077# ) and <S> E ~ tp(= 0.1903^/^) for 
Case (i), and Z E ~ -H{= -51692# ) and <5> E ~ ip(= 1.074 VZ/k) for Case (ii) at the 
starting timef^. 

The evolution of A E is shown in Fig. 1. The solid line represents Case (i), while the 
dashed line represents Case (ii). Although A E with all the constraints grows in almost 
the same way as in the Einstein de-Sitter case, that without 3) is greatly enhanced by 
a factor ~ 6 x 10 2 . To see why, consider the time evolution of gravitational "constant" 
G and its change rate G, as shown in Fig. 2 and Fig. 3, respectively. Corresponding 
to the enhancement in the growth rate of A E , we see a great deviation of G from 
G\ P resent in Case (ii). On the other hand, the variations of gravitational constant itself 
are within about several percent in both cases. This means that it is not the variation 
of gravitational constant itself but its change rate that is the key to the enhancement of 
the growth rate of perturbations. 

We emphasize that the enhancement of A E shown above is a generic feature for a wide 
range of parameters, such as Qdustfi or \G / G\ presen t. In order to see that the qualitative 
behavior is the same for the case of &dust,0 7^ 1-0, consider Fig. 4, in which Qdust,o 
(). 1 . We set \G/G\ present = 2.747 x lO^yeai-- 1 for Case (i), \G/G\ present = 2.889 x 
10" 14 year _1 for Case (ii), and set at the starting time Z E = -19138F , $ E = 1.186^3/^ 
and Z E = -10023F , ® E = 0.1228^/3/^ for each case. The enhancement is slightly 
reduced because of the lower density of matter. As for the values of \G/G\ presen t, as we 
mentioned above, constraint 3) strictly confines them. In order to satisfy constraint 3) in 
the BD theory, \G/G\ presen t must be almost exactly \G / G\ presen t = 2.039 x 10~ 14 year _1 
for ttdustfl = 1-0 and \G/G\ present = 2.747 x 10 _14 year~ 1 for U d ust,o = 0.1. The other 
values correspond to large deviations | AG /G\ nuc i eosyn thesis, and A E is greatly enhanced 
in most of such cases. 



5 Concluding Remarks 

In this article we have explicitly presented gauge-invariant cosmological perturbation 
equations in generalized Einstein theories which are characterized as equations ( |3.1| ) and 
( [A. 21 ) with a perfect fluid. Although such equations have been of interest to resolve the 



galaxy formation problem, the complex structure of the system has made the formulation 
difficult. Therefore, we have first simplified the fundamental equations through conformal 
transformation. This process reduces the equations to the same form as in the Einstein 
theory which allows us to apply many well-developed techniques to their solution. This 
considerable simplification is important not only for this problem, but for several others 
as well. Then we have applied the covariant approach to derive the perturbation equa- 



3 We include the minus sign to pick up the growing mode. 
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tions. This approach makes the formulation straightforward, and the intuitive meanings 
of the naturally introduced gauge-invariant variables are easy to obtain. Now that its 
correspondence to the conventional Bardeen's approach is well understood, this approach 
is becoming established and will be applied in various situations to formulate relativistic 
perturbation theories. 

As a simple example, we have applied the derived equations to perturbations in 
Brans-Dicke cosmology. We have presented an analytic solution for a specific background 
model, and shown an enhancement in the growth rate of perturbations. We have also 
analyzed perturbations in more generic models numerically. The enhancement of density 
perturbations occurs throughout the range of the initial conditions or parameters, so 
we have evaluated them assuming the observational constraints. Although the three 
constraints we consider here are too strict to permit a large enhancement in density 
perturbations, if we take into account only today's direct observations, we have found a 
growth rate high enough to account for galaxy formation. 

Our calculations imply that it is not the variation of effective gravitational constant 
itself but its change rate that is the key to the acceleration of the perturbation evolu- 
tion. That means that if only our universe experienced an epoch in which the change 
rate of gravitational constant was high, even if the variation itself was small, then den- 
sity perturbations could have grown fast enough to form the structure observed today, 
and that other complex processes would be wholly unnecessary. Hence the structure 
formation mechanism in generalized Einstein theories is very attractive. It may also be 
difficult to rule out completely such a possibility from observations. If GETs correctly 
describe nature, they may predict enhanced growth rates for density perturbations. We 
are now engaged in further studies of theories with non-minimally coupled scalar fields 
and perfect fluid matter, to which we plan to apply the results presented here. 
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Appendices 

A Formulation for the Most General Type of GETs 

In this appendix we consider the most general action for GETs, i.e., 




(Al) 
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where F((j), R) is an arbitrary function of a scalar field <j) and of a scalar curvature R. 
Here, we do not include Case (A) (where F is a linear function of R) because it has been 



discussed in the text. Assuming a perfect fluid as (|3.2| ), the basic equations arep6| 



G, 



ab 



+ 



\dRJ 
V a V 6 



V a( />V 6 0-- 5afe (V0) 2 



+ ^9ab ( F 



dF 
~dR 



R 



dF\_ n fOF\ 



dR 



\dR 



dab 



| 1 rpM 

+ 2 ab 



(W) 2 + 1? = o. 



(A2) 
(A3) 



Although these involve higher-derivative terms of the metric, the conformal transforma- 
tion dU) with 

(A4) 



u = - In | 2ir 
2 



<9F 



dR 



eliminates those higher-derivatives, and the introduction of a new "scalar" field 

dF 



Kip = \J~Qu) = —— In ( 2/r 



dR 



(A5) 



enables us to write down the equations in the Einstein form: 

Gab = ^ haf^b^P ~ ^ab(V^) 2 ) + K 2 E((j), if) ( V a (/>V 



where 



1 



-k 2 U(^ p)g ab + K 2 N{p) [nu™ +ph%} , 
zarf) 3 



E(<t>,cp) 
N(<p) 



(sign)e(0) exp 



-KLf 



(sign) 



-Kp 



OF 
dR 



dF 



dF 



dR 



(sign) exp 
dF 



N(<p) 



dR 



~V6kp 



R-F 



M 



M 



h 



a 

Mb 



"a ' 



h M 
n ab 



= e~"u a M , 



tab 



-2uiiab 

e n M . 



(A6) 
(A7) 

(A8) 

(A9) 

(A10) 

(All) 
(A12) 



The (sign) equals 1 for the present universe, and for a generic spacetime there appears 
a singularity at dF /dR = beyond which the (sign) turns to be —1. Then we should 
usually assume (sign) = 1 although we leave it in this article. The equation for ip is 
obtained by taking the trace of equation ( A2 ) and using relation ( |A5[ ): 



6 



dF 



Up + ^rKE{4>, p)(V<P) z + ^rKN(p) \ —R -2F --(-» + 3p) } = 0. (A13) 
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In these equations, R and F are functions of 4> and <p i.e., R = R((f>, <p) and F = F(cj), ip) 
F((f>, R((f>, <p)), which are defined implicitly by the relation 



dF 



a given function of <f> and R 



(sign) ^ exp 
2k z 



(A14) 



Tab 
p* 



V* u a u b + P* h ab + 2*(a«&j + ^*ab, 

hp 2 + ^(0, ip)^ 2 + y>) + N(<p)n, 

Itf + 1^(0, ^ _ Ufa ^ + N ^ pj 



o 



The gravitational field equation ( |A.q ) corresponds to the Einstein one with two scalar 
fields and a perfect fluid interacting with each other through E(<fi, ip), U(4>, (p) and N((p). 
The energy- momentum tensor in equation ( [A. 6 ) is expressed in a fluid form like equation 
([DJD as 

(A15) 
(A16) 

(A17) 
(A18) 

K*ab = 0, (A19) 

(A.20) 



v; 



l M 



l a i 



to first order. Hence, the derivation of the GI perturbation equations to follow is straight- 
forward, and proceeds as described in §3 

Particularly in the case of F = L(R) and e(4>) = (Case (B)), namely, when the 
action is given by 

S = J d 4 x^9 \L(R) + L m ] , (A21) 



the basic equations become 

G n h = 



k 2 \ V a ipV b <p - -5afe(V(^) 2 - U(ip)g ab 



+ 



K 2 N{<p){nufu¥+phH} 



□ <p 



dU(ip) 
dip 



+ n 2 M((p){fi - 3p} = 0, 



with 



N(<p) 
M(ip) 



(sign) exp 
dL 



N{ip) 



dR 



-VSnip 



R-L 



(sign) 



dL 

dR 



dL 



dR 



OK 



(A22) 
(A23) 

(A24) 
(A25) 
(A26) 



Here, R = R(<p) and L = L(ip) = L(R(ip)) through the relation 

dL 



dR 



(R) 



a given function of R 



(sign) exp 



V6 



Kip 



(A.27) 
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Ill 



Equations flA.22 ) and ( A.23|) have exactly the same form as equations (|3.7| ) and (]3[ 
§|3[ Hence the perturbation equations are given by equations ( B) , (|3~50l) and (^49| ). 
Among the models in Case (B), we are mostly interested in R 2 gravity from the 



astrophysical point of view. By setting L(R) 



( |A.24|) ~ (|A.27D become 



with (sign) = 1, equations 

(A.28) 
(A.29) 

2 

(A.30) 
(A31) 



In |28|| , perturbations in the R 2 gravity were discussed. However, they were concerned 
only with metric perturbations without matter fluid (L m = 0), and their approach was 
not GI, although their gauge (conformal-Newtonian gauge) gave physical modes. 



N(<p) 


= exp -^VGklp , 




M(<p) 


= ^AV). 

OK 




C%) 


= Sk^Y 6XP 


~Vg 


Kif 


= ^ln[|l + 2Cfl|]. 



B Perturbations in a Scalar Field Dominated Universe 

In this appendix, taking the limit fi, p — ► 0, we show our perturbation equations reduce 
to the ones in a scalar field dominated universe derived by Bruni et al.[22]. Since they 
assumed the Einstein theory, we must set 



1 

2k 2 ' 
h 



in the action fl3.1| ), which immediately gives 

N{<p) = 1, 

A%) = 0, 

[%) = V{4>). 



(B.l) 
(B.2) 

(5.3) 
(5.4) 
(5.5) 



with ip = <p. Besides, as they used center of mass frame u E , we can use equations 
0)~0) as the basic equations. 

In a scalar field dominated universe, matter fluid perturbations are absent, so we 
need consider only two independent GI perturbations. Bruni et al. adopted the total 
energy density perturbation A E = (a 2 / fx^^'V 2 ^^* and the curvature perturbation C E = 



a 
as 



2(3) V 2 



^R 



for them. They are expressed in terms of our GI perturbation variables 



A 



C 1 



-4a 2 HZ E + 2K 2 a 2 n*Af. 



(5.6) 
(5.7) 
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We can recover their basic equations in the limit of /i, p — ► in our equations. The 
background equations ( 3.31 )^( ^.34 ) become 



35 + 35 2 + k 2 (c/> 2 - V(<p)) = 0, 
or 

+ 3i?0 + y'(0) = o, 

where /U* = i(/> 2 + V{4>) and equation ( . 1 0| ) is equivalent to 

/i* = — 350 2 . 



(B.8) 
(B.9) 
(5.10) 

(5.11) 



The condition that total energy flux vanishes in the center of mass frame implies in this 
limit 



(B.12) 



which gives <5 E — > 0, so that A E has a finite value even in the limit of fj, — ► 0. 

In this limit, perturbation equation ( 3.48| ) is trivial. Equations ( |3.49| ) and ( 3.50 ) give 
respectively 



2</> 2 \ a 2 a 2 



A 



3/7 + 



4V 



+ 



2V' 



/2 



+ v" 



^z 1 



a* 4> 



(5.13) 
(5.14) 



By using the relation ( B.6 ) and its derivatives, we can eliminate § E / n, <& E / [i and & E / ^ 
from equation ([B.14 ) as 



A 



-35 + 



350 2 K 



2x2 



2H 



Af + 



4a 2 5/x* 



C 1 



(5.15) 



Rewriting 35 for and using the notation 7 = 2 / '/i*, this becomes 



A 



37 
4a 2 # 



(7 - 1)0 



3K 2 /i*7 
20 



(5.16) 



which corresponds to equation (55) in [^| . Moreover we can eliminate Z E from equation 
( |5.13| ) by using equation ( 5.7|) as 



2K 2 a 2 /i*Af 



8K 2 a 2 Hfi^ + 



V)^ 
2ac , "2 



(0 2 -V)- 6kV5<^ + 2(-2n 7 + 6k) 



35 



5/i* 



Af. 



(5.17) 
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Replacing by equation ( |B.16j ) and by using equation (|3.47| ), this becomes 

C E = %C* + l^(3)v 2 Af + 4k( 9 — *4f) Af , (5.18) 



which corresponds to (56) in Ifffil . The equations (|£>.16|) and (|£>.18|) are the basic equa- 
tions for GI perturbation variables A E and C E given in [22]. 



C Analytic Solution in a BD Cosmology 



The solution (|4.2| )^ (|4.4| ) of the flat dust BD universe with the condition (|4.l| ) is described 
in terms of conformal variables as 



t 



5+3^ 

4 + 3ul ( t\ 4+3u 



5 + 3uj \to 
V2(3 + 2u) 
k(5 + 3a;) 



In 



5 + 3a; t 



4 + 3to t 



5 + 3a; i \ 5+3u 



4 + 3uj t 



W« 3 { — 



6(l+ui) 

2(3 + 2u) fb + 3u t \ 5+3w 



n 2 t 2 (A + 3oj) l4 + 3wt 



(CI) 
(C2) 

(C3) 

(C4) 



Inserting this explicit background solution into equations (4.15)~ (4.17), we present the 
perturbation equations as 



A 



-Z E + 



3(7 + 5lo)k 



V2(3 + 2w)(4 + 3uj) 



4> E + 



Kn 2 to 



5 + 3uj t 



y/2(3 + 2u) V 4 + 3wi 



(C5) 



(3 + 2o;)(4 + 3a;) 1 g 

(5 + 3a;) 2 t 2 
2(3 + 2a;) 1 E 
(5 + 3oj) t 

V2(3 + 2w) k 

2(4 + 3a;) I ~ ^2(3 + 2a;) 
4 + 3a; 



Kn 2 to 



5 + 3a; t 
4 + 3a;^ 



+ 2J2(3 + 2w) 



(5 + 3a;) 2 £ 2 



(C6) 



2(4 + 3a;) l^ E 
5 + 3a; £ 

4(4 + 3a;) 2 1 



2 4 + 3a; / 5 + 3a; t 
(5 + 3a;) 3 t 2 ' 5 + 3a; I 4 + 3a; ^ 



2(3+2^) 
5+3u 



+ n 
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^2(3 + 2w) 



(5 + 3a;) 3 K P 
(4 + 3a;) 1 zE 



(5 + 3u) 2 d 

Then we can solve these equations analytically as 

2(2+u;) 



A £ = ci < 



t \ 5+3t " 4(8 + 5w) /5 + 3aA 5+3^ 1 
W + (5 + 3^)2 U + 3aJ t 2 n 2 



/ 5 + 3u; A 5+3w 
r!2 \4 + 3a;^J 

4+3l^ 4+3^ 
/ £ \ 2(5+3w) / £ \ 2(5+3w) 



_ E (3 + 2w)(4 + 3w) 1 / t \ 5+3t " 

-Z = — ci- ' 



+ C2to 



(5 + 3^)2 t \t J 
( 5jh3a; t_\ 

10+7uj 



2+u 



1+ 6+5^ 
2 +u , • • 10 + 7 - 



. n /5 + 3uAB+ssr / M 2 ( 5 + 3 "> 
+ =4 7 7^ 7- J-u-iinT) 



A\4 + 3ujJ U 



1. 1 .,■ - - - — 



V v ^ + q,, I 4 , q. . ; Kn 2 



(5 + 3^)2 \A + 3uJ kuH\ 

4+3w 4+3^) 



where Ci, C2, c 3 and C4 are arbitrary constants, and 

v _ 4 + 3w 4 + 3w / 5 + 3a; t \ 

" = 2(2 2 + lo U + 3o;tW 
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Figure Captions 

Fig. 1: 

Density perturbations in the BD theory for the case of Qdusto = 10. The solid line 
represents Case (i); all the observational constraints 1), 2) and 3) are satisfied, while 
the dashed line corresponds to Case (ii); only 1) and 2) are satisfied (see text). The 
density perturbation in the Einstein theory with the same initial condition coincides with 

CO- 



Fig. 2: 

Evolution of effective gravitational constant (G=l /</>). Each line corresponds to the line 
in Fig. 1. bearing the same number. 



Fig. 3: 

Evolution of change rate of effective gravitational constant. Each line corresponds to the 
line in Fig. 1. bearing the same number, and Hq = lOOkm/sec/Mpc. 



Fig. 4: 

Density perturbations in the BD theory for the case of Qdustfl = 0.1. The solid line 
represents Case (i); all the observational constraints 1), 2) and 3) are satisfied, while 
the dashed line corresponds to Case (ii); only 1) and 2) are satisfied (see text). The 
density perturbation in the Einstein theory with the same initial condition coincides with 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9404023vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9404023vl 



